March 4, 2013 



Scaling Properties of the Lorenz System 
and Dissipative Nambu Mechanics^ 

Minos Axenides 1 and Emmanuel Floratos 1 2 

1 Institute of Nuclear and Particle Physics, N.C.S.R. Demokritos, 
GR-15310, Agia Paraskevi, Attiki, Greece 
2 Department of Physics, Univ. of Athens, GR-15771 Athens, Greece 

axenides@inp.demokritos.gr; mflorato@phys.uoa.gr 
Abstract 

In the framework of Nambu Mechanics, we have recently argued that Non-Hamiltonian 
Chaotic Flows in i? 3 , are dissipation induced deformations, of integrable volume preserving 
flows, specified by pairs of Intersecting Surfaces in R 5 . In the present work we focus our 
attention to the Lorenz system with a linear dissipative sector in its phase space dynamics. 
In this case the Intersecting Surfaces are Quadratic. We parametrize its dissipation strength 
through a continuous control parameter e, acting homogeneously over the whole 3-dim. phase 
space. In the extended e-Lorenz system we find a scaling relation between the dissipation 
strength e and Reynolds number parameter r . It results from the scale covariance, we 
impose on the Lorenz equations under arbitrary rescalings of all its dynamical coordinates. 
Its integrable limit, (e = 0, fixed r), which is described in terms of intersecting Quadratic 
Nambu "Hamiltonians" Surfaces, gets mapped on the infinite value limit of the Reynolds 
number parameter ( r — > oo, e = 1). In effect weak dissipation, through small e values, 
generates and controls the well explored Route to Chaos in the large r- value regime. The 
non-dissipative e = integrable limit is therefore the gateway to Chaos for the Lorenz system. 



^The present work is dedicated to the memory of Prof. J.S.Nicolis 
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1 Introduction 

In the framework of Nambu mechanics we have recently [lj [2] proposed a geometric approach 
for the study of dynamical dissipative systems [3J. It is implemented through a special 
decomposition of the flow of the system in two parts. The non-dissipative and threfore phase- 
space volume preserving component and the dissipative constituent part which is volume 
contracting. In the particular case of D = 3 phase-space dimensions such a splitting takes 
the suggestive form 

'x = ViTi x VH 2 + VD (1.1) 

where H\ , H2 are the Nambu "Hamiltonians" or equivalently, the Clebsch-Monge potentials 
[I] where D is the dissipation potential. In particular, for the Lorenz system [HE], 

x = a(y — x) 

y = x (r- z )-y (1.2) 

i = xy — bz 

the functions H\,H.2,~D are: 

Ht = \(y 2 + {z-rf) 

H 2 = az-^ (1.3) 

D = - - (ax 2 + y 2 + bz 2 ) 

The non-dissipative Lorenz system which was studied carefully in pQ was also investigated 
some time ago in [7] and in a similar vein further back in time by [8] 
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'x = ViTi x VH 2 (1.4) 

is an integrable dynamical system, where the trajectories in phase space are given by the 
intersection of the surfaces defined by the conserved Nambu " Hamiltonians" (H\,H2) . In 
this picture the dynamical system in rel(l.l) is defined through its non-dissipative part of 
eq,(1.4), over which its phase space dissipation V-D operates. In order to control its dissipation 
strength and take hold as accurately as possible of its phase space volume contracting effect, 
we introduce an associated control parameter e as follows: 

f = V#i x VH 2 + eVD (1.5) 

It is defined in a range of values on the interval [0, 1] with its value e = recovering the Lorenz 
system nondissipative part, whereas for e = 1 the full system being obtained. Clearly the 
e-parameter controls the dissipation rate of the phase-space volume of the dynamical system 
in question. For the Lorenz case, for example, 



SV(t) = e -< 1+a+b ^ 6V(0) (1.6) 

We could say that this is nothing else but a rescaling of time. As we shall demonstrate in 
the particular but very interesting case of the Lorenz system, the rescaling of time is possible 
only under appropriate simultaneous rescalings of all the dynamical variables x,y,z of the 
phase space, the time t and more importantly, of the Reynolds number r , while keeping the 
two others b, a intact. 

At this stage we could observe that the introduction of the control parameter e measures 
also the strength of the attraction of the Lorenz ellipsoid [5]. 

Indeed, the rate of change of the Liapunov function 

, --(rx 2 ' ~ 2 ' ~'~ *-*\--™. i ~- ■ 3 ~ 2 - 
2 V 

is given by 



V = - [rx 2 + ay 2 + a(z - 2r) 2 j = oH x + rH 2 + -r 2 a (1.7) 



V = r- W = -&T 



rx + y + b(z — r) — br (1. 



When e is a small positive number, the velocity is entering the Lorenz ellipsoid V, if all 
of its points are located outside the ellipsoid defined by the vanishing of the expression of the 
last equation. On the other hand, in the case of the nondissipative system of eq.(1.4), i.e. 
e = the velocity of the particle is tangent to the ellipsoid V = V(xo, Uoi z o) where (xq, yo,zo) 
are the initial conditions of the motion, because H\,H 2 are conserved. Thus e controls the 
overall attractive strength of the attractor. 

In our recent work we also introduced a method of Matrix- Heisenberg Quantization for the 
Lorenz system aiming to examine the compatibility of the classical Lorenz strange attractor 
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dynamics with the fundamental principles of Quantum Mechanics. We defined rigorously 
the quantization procedure starting from the non-dissipative system (1.4) where we replaced 
the phase space coordinates x,y,z with N x N hermitean matrices X, Y, Z consistently with 
the appropriate rules of Quantum Correspondence. The quantum dissipation was defined 
to respect the above quantum correspondence principle. It is intuitively obvious that, if 
dissipation strength in phase space is big enough, the approach of the system to the classical 
limit is fast and efficient expressed through the exponentially fast vanishing in time of the 
commutators of X, Y, Z. Thus it is imperative for the study of the coexistence of strange 
attractors and quantum mechanics long time scales to be maintained in the non- vanishing of 
quantum commutators. This is feasible through a continuous dissipation strength controlling 
parameter e. 



2 Controlling Dissipation in the Lorenz System 

Before we study in detail the scaling properties of the extended Lorenz system ( we shall 
refer to it as the e-Lorenz system) 



x = o(y — ex) 

y = x {r- z )-ey (2.1) 
z = xy — ebz 

lets analyze the stability properties of its critical points. On this issue we keep in mind 
that a similar model has been considered in the literature ( [911101 [12]) but with precisely fixed 
control parameter e = ^7=, in order that the limit r — > 00 becomes explicit. For the stability 
analysis of the critical points we follow [6] and [12] . 

The critical points of the system of eq. (2.1) are: Po : (x = y = z = 0) the origin of the 
phase space, and P± : (x± = ±y/b(r — e 2 ), y± = ±ey/b(r — e 2 ), z± = r — e 2 ). In order that 
P± are real points, we must have r > e 2 . Linearization of the system around Pq gives 
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provides the characteristic Polynomial 



fo(A) 



-{X + be) [(ea + \)(e + \)-ar\ 



(2.4) 



It posesses two negative and one positive eigenvalues. They correspond in turn, to a 
stable 2-dim. manifold Wq with two attracting directions in i? 3 as well as to an unstable 
1-dim. manifold Wq with a single repulsive direction, Vu > 0, r > e 2 . 



Ai 
A 3 



-be 



e{l + a) + Je 2 {l + a) 2 + \o(r - e 2 



e 2 (l + <r) 2 + 4cj(r-e 2 ) - e(l + a) 



(2.5) 



Linearization, around the P± critical points, gives identical characteristic polynomials 



p ± (A) = A 3 + e(l + a + 6)A 2 + b{r + ae 2 )X + 2a6e (r - e 2 



(2.6) 



Since all the coefficients are positive this cubic polynomial has at least one negative 
eigenvalue(attractive direction in i? 3 ). 

The bifurcation from 3-dim stable manifold (dim W± = 3) to 1-dim stable manifold ( 
dimW^ = 1 ) occurs, when the real part of the pair of complex conjugate roots A = Xr ±iA/ 
crosses the value Xr = 0. Simple algebra shows that Xr = is equivalent with the statement 
that the constant term in P± in rel.(2.6) is the product of the coefficients of the linear and 
the quadratic terms. Thus the critical value of r is given by, 



r c (e) 



e 2 a 



a + b + 3 



(2.7) 



a-b-1 

For r > r c (e) the subcritical Hoph bifurcation leads to the aperiodic attractive orbits of 
the famous Lorenz Strange attractor[5]. From ref.(2.7) we get the interesting result that for 
small e close to the non-dissipative integrable limit of e = the strange attractors appear for 
small values of r > r c (e). 

In concluding we present for the critical case of Xr = 0, the other roots of the cubic 
polynomial P±. The real root is given by 



A = -e(l + a + b) (2.8) 
whereas the imaginary component of the complex conjugate one is: 

_2eV(a + l) 
A/ " a-b-1 (2 - 9) 

Xj represents the angular frequency of the orbit rotation around the critical points P±. 
In what follows we are going to investigate the scaling properties of the extended e-Lorenz 
system which will make transparent the e-dependence of r c (e) in rel. (2.7). 
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3 Scaling Relations in the Lorenz Model 

Firstly we rescale independently all the variables of the system a, /3,7, A G R 

x = ax' , y = j3y , z = 72/ , t = Xt' (3-1) 

We demand that in the transformed primed system a similar structure of constants appear 
in the equations of motion, i.e. of the r, e, b, a type ( the dot represents time derivative with 
respect to t'), 

x> = a '(y'-e'x') 

V' = x'(r'-z')-e'y' (3.2) 
z = x y — e z 

we obtain 

e' = Ae, b' = b, a' = a, r' = X 2 r, a=\, ^ = 7 =^ (3-3) 

A A 

The dynamical variables x, y, z are related to their primed ones as: 

x(t) = j x'(Xt) , y(t) = ^y'(Xt) , z(t) = ^z'(Xt) (3.4) 

We will show below that the established covariance property of the e-Lorenz system under 
the scalings (3.3 , 3.4 ) implies interesting constraints for the dependence of the phase space 
coordinates x, y, z on the time, the parameters r, e, a,b as well as on the initial conditions 

Indeed the time evolution of the phase space coordinates is given by the exponential of 
the Liouville operator Co, (xo = x(0),yo = y(0),zo = z(0)), 

x(t;r,e,a,b,xo,yo, z ) = e tCo ■ xo (3.5) 

where 

C = a (y a - ex ) d XQ + [x (r - z ) - ey ] d yo + (x yo - tbz ) d zo (3.6) 
Under the scaling (3.3-3.4) which we rewrite explicitly as 

t ->■ , cr ->■ a , b^b , r ->■ X 2 r , e ->■ A e (3.7) 
A 

as well as, 

xo ->■ Axo , yo ->■ A 2 yo , z ->■ A 2 zo (3.8) 
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( VA G R ), it is straightforward to check that Cq scales as 



£ -> A£ (3.9) 
At last from rel.(3.5) we obtain that Vt G i?, 



A 2 r, Ae, a, b, Xx a , X 2 y a , \ 2 z Q ) = A x(t; r, e, a, b, x , y a , z a ) (3.10) 



and for the y(t) and z(t) coordinates we obtain respectively 

y -»• A 2 y , z^A 2 z (3.11) 

We note that the nondissipative Lorenz system (e = 0) satisfies scaling relation (3.10). 
When (e )we impose A • e = 1 on the LHS of eq.(3.10) which implies that the whole r - e 
plane is foliated by the continuous set of parabolas 

T 

~2 = rc = constant (3.12) 

Each parabola corresponds to a different phase of the Lorenz system which is specified 
by the value of the Reynold's parameter r for which the parbola cuts the e = 1 line . On the 
other hand if we fix a given parabola all of its points are related through the scaling relations 
(3.3-3.4). Indeed the condition A • e = 1 implies that: 

( r 1 1 1 \ 

x(t; r, e, a, b, x , yo,z ) = e x [ e t; ^, 1, a, b, -x , -^yo, ^zq) (3.13) 



Similarly for y, z we have 



i r 111 

y(t; r, e, a, b, x , y ,z ) = e 2 y ( e t; 1, a, b, -x , -^yo, -^z ) (3.14) 



( r 1 1 1 \ 

z(t; r, e, a, b, x , y ,z ) = e 2 z f e t; ^, 1, a, b, -x , -^yo, ^oj (3.15) 

Equations (3.13-3.15) constitute the main result of the present work. It relates the time 
evolution of the extended e-Lorenz system e / 1 with the standard Lorenz one e = 1. Indeed, 
by denoting with rc the ratio rc = ^ we can read off eqs.(3.13-3.15) in two distinct but 
complementary ways. 

Firstly, if we fix rc we see that the r parameter of the e-Lorenz system ( LHS of eqs.(3.13- 
3.15)) ranges between the values of < r < rc when < e < 1. It follows that the physics 
of all the points on the parabolas r = rc ■ e 2 of the ( r - e ) plane is the same. 
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Chaotic Lorenz Zones 




Period Doubting Bifurcations 
Intermittent 



Homoclinic Explosion Orbits 
Lorenz Strange Attr actors 
Preturbulence 



0.0 0.2 0.4 6.6 0.S 1.0 



Figure 1: Lorenz System Structure of Dynamics in r-e plane 
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On the other hand if we fix r for the e-Lorenz system ( LHS of (3.13-3.15) ), by varying e 
(0 < e < 1), we are able to scan the region of the standard Lorenz system, for rc, r < rc < oo 
(see fig.l). 

Some remarks are in oder: 

By fixing r in the e-Lorenz system and taking the limit e — > we recover the non- 
dissipative Lorenz system of eq.(1.4). It is integrable and describes the motion of a particle 
in an one dimensional anharmonic potential jl]( or the pendulum). We can still rescale this 
sytem by choosing A = ^ and we discover the infinite r limit of the full Lorenz system, which 
has been studied in detail in ref. [10|. [TTj . In this limit the Lorenz attractor is degenerate into 
an 8-figure stable limit cycle ( see also in |14j). In the above mentioned works, the - correction 
has been studied and simple bifurcations of the limit cycle have been observed. In [10] an 
exhaustive search has been made by lowering the values of r. It has been discovered that 
for relative small values of the 8-figure still survives somewhat deformed nevertheless up to 
r = 300. 

In table 1 we summarize the behavior of the Lorenz system for various interval values of the 
Reynold's number rc which comprise the famous Feigenbaum Route to Chaos. In the present 
case of the extended e-Lorenz system with a (r, e ) space of parameters they correspond to 
parabolic zones of equivalent dynamical systems, specified by the scaling relation r = r^e 2 . 
In table 1 we depict two such characteristic equivalent Lorenz systems for e = 1 and e = .01. 
The Route to Chaos for each of them is characterised by the following distinct behaviours 
for different values of r for fixed e : 

For values of the r = rc with e = 1 in the range (30 < r < 313) is composed of an 
infinite cascade of period doubling bifurcations (2, 2 2 , 2 3 , • • • , 2 ra ) which appear until a critical 
value of r , r c = 30.1 is reached below which the strange attractors appear for r in the 
range 24.06 < r < 30.1 . At approximately the value of r = rc = 166.07 we have the 
Intermittency scenario to Chaos which was first described by Pomeau and Maneville |15j . 
For 13 < r < 24.06 an interesting "preturbulence" regime has been observed in the works 
of J.Yorke et.al.|16j. For a complete and detailed description of the behaviour of the Lorenz 
system for all values of r one should look at the work of Sparrows who docs it exhaustively |6j. 
All of the above are depicted in figure 1. 
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Table 1: (r-e) Lorenz Roadmap to Chaos 





r = r c , e = 1 


r , e = 0.01 


Period Doubling Bifurcations 


59.5 < r < 313 


0.00596 < r < 0.0313 


Intermittency 


r = 166.07 


r = 0.0166 


Orbits from Homoclinic Explosions 


30.1 < r < 59.5 


0.0030 < r < 0.00596 


Strange Attractors 


24.06 < r < 30.1 


0.0024 < r < 0.00301 


Preturbulence 


13.9 < r < 24.06 


0.00139 < r < 0.0024 



4 On the Scale Invariant Lorenz System 

In the last part of this work we introduce a form of the e-Lorenz system which is by construc- 
tion invariant under the scalings in rel. (3.3-3.4). To this end we define new independent and 
dependent variables ( r > e 2 ) : 



Vr - eH , X 



x 



r — e 



2 ' 



Y 



r — 



r — 



(4.1) 



which satisfy the system of eqns. 

X = a(Y - (X) , Y = X(l + C 2 - Z) - (Y , ( = XY - b(Z (4.2) 
with C = / — — — *t 

The derivative "dot" is with respect to the new time r . Every parabola ^ = c in the r - 
e plane is determined by a fixed value of 

1 



or 



c 



r = (1 + ^2 |f" 



1 



(4.3) 



(4.4) 



By construction it is easy to see that the variables r, X,Y,Z are scale invariant under rel. (3.3- 
3.4). We may also observe that in the system (4.2) the parameter ( appears in both places 
where the Reynolds number as well as the dissipation control parameter eappeared in eq.(2.1) 
independently. 

Repeating the critical point analysis, as before, and their stability we find that there are 
three critical points, the point P D = (0, 0, 0) with two attracting and one repelling directions 
with corresponding eigenvalues: 



Ai = -6C, A 2 = --[C(l + <r) + ^/C 2 (l + <r 2 + 4<7] , 
A 3 = \ K/C 2 (l + ^ + 4a) - C(l + *)] 



(4.5) 
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as well as two symmetric , under the reflection (x, y, z) 
P± = (±Vb,±(Vb , 1) 

with characteristic polynomial 



-x, —y, z)critical points 



P±(X) = X 3 + ((l + b + a)X 2 + b l + ( 2 (l + a) A + 2b(a 



(4.6) 



Once more the bifurcation condition for the appearance of the strange attractor is £ < Q c 
where 



c; 2 



(7-6-1 



(<T + l)(ff + 6+l) 



For the standard values of a = 10 and b = | we find 



Co 2 



19 11 



(4.7) 



(4.8) 



11-41 451 

We note finally that the scale invariant form of the evolution eqs. (4.2) can be cast in a form, 
which exhibits the role of dissipation on the non-dissipative sector of the system. In figure 
2 we plot rc = \ ,i- e - the relation of rc against Each value of Q determines a unique 
parabola in the r-e parameter plane of the equivalent physics of the Lorenz system. 




Figure 2: Lorenz ( Dissipation Parameter 

The Non-Dissipative part was defined in rel.(1.4) and corresponds to the £ = case : 

X = aY, Y = X(1-Z), Z = XY (4.9) 

We choose to work with the Cylinder intersecting with a Paraboloid as the two Nambu 
"Hamiltonians" among the whole set of SL(2,R) geometries pQ. 

X 2 



Hi 



1 



y z + (i - zy 



H 2 = aZ 



or 



X = Vifi x VH 2 



(4.10) 



(4.11) 
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For ( = the two surfaces Hi , H 2 are conserved 



Hi =fli(t = 0) , H 2 = H 2 (t = 0) 



(4.12) 



and they define by their intersection the trajectory of the system (4.9). When £ 7^ 0, Hi,H 2 
are not conserved. Still one can use in the place of X, Y, Z the variables X and H 2 in order 
to elucidate what happens. For (/Owe obtain [TT] 

"X 2 



X + (l + a)£X + X 



+ H 2 



a 







and 



#2 + c 



bHo 



2a' 



)X> 







(4.13) 



(4.14) 



We see now the crucial role of ( in the introduction of two terms: Firstly a friction term 
proportional to X and secondly, a memory term in the anharmonic potential minimum By 
eliminating H 2 and using rel.(4.14) we obtain the Takeyama evolution memory term, which 
changes in a non-Marcovian random manner the symmetry of the single well potential into 
the double well one [TT] . 
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5 Conclusions-Open Problems 



The introduction of the e-Lorenz system, which controls the strength of the dissipative sector 
implies the existence of a weak (under-dissipated) and strong(over-dissipated) phases for the 
Lorenz system. They are separated in the r — e parameter space phase diagram by the critical 
line e = 1 which reproduces for different values of r the route to chaos of the Lorenz system, 
which is presented so lucidly in the work of Sparrow[6j. We have demonstrated therein that 
a Lorenz system with weak dissipation ( small £ values) is equivalent to the one with large 
values of the Reynolds number r for the standard Lorenz system with e = 1. 

The exploration of the scaling properties though, brought unexpected new information 
about the e Lorenz which smoothly joins to the e = integrable case. The scaling eqs.(3.3- 
3.4) show that we have very specific combinations of the time, the initial conditions and the 
parameters r and e in each order in the Taylor expansion of the solution with respect to time, 
such that the scaling relation holds true. Also the bifurcation behaviour of the solution can be 
controled by the one parameter Q = ~^^g in rel.(4.3) which parametrizes the foliation of the 
r-e plane by one parameter family of parabolas. All the points of each parabola are physically 
equivalent. Last but not least, the Nambu surfaces being the appropriate geometrical tool 
for the C = integrable case, are useful also for the small £ range which corresponds to the 
large r model (e = 1) where we know that successive bifurcations of the figure 8 periodic limit 
cycle |10t ITT] lead to the strange attractor configurations. The Non-dissipative Limit is thus 
the gateway to Chaotic and turbulent Flows for the Lorenz system. Interestingly, it may be 
also the gateway to Quantum Chaos in a matrix formulation of the e-Lorenz system, where 
the expected presence of decoherence can become suppressed in an controllable way. This 
interesting possibility will be investigated in a separate work in the near future. 
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